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Abstract 

In this paper, we calculate the transverse momentum diffusion coefficient of 
heavy quark and gluon penetration length in the deconfinement phase of Sakai-Sugimoto 
model, which is known as a holographic dual of large N c QCD. We find that for the 
heavy quark moving through the thermal plasma with a constant velocity v < 1, the 

transverse momentum diffusion coefficient kt oc A73 J 4 / Td, and the gluon penetration 

2 

length Ax oc E 5 . These results are different from those calculated in N = 4 super- 

1 , .i 

Yang- Mills theory, which are kt oc A72 and Ax oc E , respectively. In the high 
energy limit, the difference between the two pairs of results should be evident, so we 
hope that the future LHC experiments can tell us which model is more closely related 
to the realistic strongly coupled QCD at finite temperature. 
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1 Motivation 



The experimental relativistic heavy ion collisions have produced much evidence signalling 
that Quark Gluon Plasma (QGP) has been formed at the Relativistic Heavy Ion Collision 
(RHIC) [U [2J. One piece of strong evidence is that the production of the high-p r particles 
is suppressed [31 H]. To explain this phenomenon within the framework of QCD is diffi- 
cult, because recently, researchers have found that QGP is a strongly coupled fluid. In the 
framework of AdS/CFT [5], one can solve problems in strongly coupled gauge theories by 
considering the corresponding problems in dual weak coupled gravity theories. So, many 
people try to solve these problems in QGP, by transferring them into a gauge theory which 
has a gravity dual and can mimic QCD to some extent. Along this way, H. Liu, K. Rajagopal 
and U. Wiedemann define the jet quenching parameter q, via a light-like Wilson loop [5]. q 
is the transverse momentum squared transferred from medium to either the initial parton 
or the radiated gluon, it is related to the average medium-induced parton energy loss by 
BDMPS formalism [24]. Meanwhile, J. Casalderrey-Solana, D. Teaney and S. Gubser prefer 
to use the transverse momentum coefficient k,t [El El HQ] . Although the two groups adopt 
different parameters, the calculations are both carried out in the same background, which 
is the y4rfS , 5 -Schwarzschild space-time. Subsequent work [TJl [131 E] includes computing q 
in backgrounds with non-zero chemical potential. Above all, the background metrics they 
use usually involve an asymptotically AdS^ component, since the gravity theory in AdS§ 
is dual to the Af = 4 SYM, their results actually apply to Af = 4 SYM. However Af = 4 
SYM is not the same as QCD, so it is problematic whether or not their results really capture 
some features of QCD, if it does, then these features should also appear in other models 
approximating QCD and having gauge/string duality, since all these models belong to one 
framework. 

Fortunately, in paper [15], Sakai and Sugimoto provide us such an new model, which we 
call S-S model in this paper. This model is a holographic dual of four-dimensional, large N c 
QCD in the low energy regime. In the high energy regime, the gauge theory in S-S model 
shows some differences from QCD, such as K-K modes. A lot of papers [16] have been 
done on this model, in which they recover some features similar to realistic QCD, such as 
confinement-deconfinement phase transition and chiral symmetry breaking-restoration phase 



2 



transitiono The calculations about screening length and jet quenching parameter q in this 
model [I2J [T5] have been carried out. But the transverse momentum diffusion coefficient 
kt has not been obtained. To give a complete comparison between above two models, we 
calculate kt in S-S model. During the preparation of this paper, S. Gubser et al pT| put 
forward a new approach to estimate the jet quenching parameter q by considering the gluon 
energy loss in the thermal plasma of strongly coupled Af = 4 super- Yang-Mills theory. Using 
this new method, we estimate q in S-S model. If we did not try this new way in S-S model, 
the comparison between the two models is still incomplete. 

This paper is organized as follows. In Section 2 we give a brief review of S-S model. In 
Section 3, after a short review of momentum diffusion constant in M = 4 SYM, we calculate 
the same transport coefficient in S-S model. In Section 4, we compute the lower and upper 
bound of gluon penetration length in S-S model and prepare to estimate q. In Section 5, we 
use the results of the previous two sections to perform a quantitative analysis. 



2 A brief review of S-S model 



In [15] . Sakai and Sugimoto present a holographic dual of four-dimensional, large N c QCD. 
This model is constructed by placing Nf probe D8-D8 into N c D4 brane background(iV/ <C 
iV c ), where supersymmetry is completely broken by compactifying the N c D4 branes on 
a circle of radius R with anti-periodic boundary conditions for fermions [20]. At low 
energy, the D4/D8/-D8 system yields a U(N C ) gauge theory with fermions, and there is 
also a U(Nj)l x U(Nf)n chiral symmetry. Being well studied in [16], this model contains 
confinement-deconfinement phase transition, the critical temperature is Td = 1/2ttR. 

When the system arrives at a temperature T < T d , the dual gauge theory of S-S model 
is in the confined phase, we should use the following background to describe it 



/ \ 3/2 / /? \ 3 / 2 

ds 2 = [-dt 2 + 5 l:j dx i dx j + f(u)dxl] + ( — j 

\RdaJ V u J 



du 2 2 , 



F( 4 ) = * = g.(-£-) , Rl^^g s N c ll f(u) = l-(^ 3 



3/4 

4 \RdaJ \« 

where t is the time direction and x l (i = 1, 2, 3) are the uncompactified world-volume coor- 
dinates of the D4 branes, £4 is a compactified direction of the D4-brane world- volume which 



1 The low as well as high spin mesons and their motions through QGP were studied in [17] . 
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is transverse to the probe D8 branes, the volume of the unit four-sphere Q4 is denoted by V4 
and the corresponding volume form by €4, l s is the string length and finally g s is a parameter 
related to the string coupling. The submanifold of the background spanned by £4 and u has 
the topology of a cigar. The tip of the cigar is non-singular if and only if the periodicity of 
X4 is 



4vr ( R m 



3 \ 1/2 



Sxa = — -™ = 2tcR. (2) 
3 \ u\ J 

When T > T d , deconfinement happens, we should use another background to depict the 
dual gauge theory, 



ds 2 = ' \-f(u)dt 2 + Sijdx'da? + dx 2 A + ' 

\Rd4/ \ U J 

O AT / \ 



u 2 dQl + — - 
f( u ) 



(3) 



This background involves a black hole. The Euclidean time direction now shrinks to zero 
size at the minimal value of u, u — ut- In order to avoid singularity, the Euclidean time 
direction must have a period of 



4vr ( R? D ^ 1/2 



3 V U T J 

QGP is the deconfined phase of QCD, in the following, we focus on the deconfined phase. 

In the deconfined phase, there exist two kinds of configurations of the probe D8 and D8 
branes. The one depicted in Figure (la) signals the breaking of chiral symmetry, and the 
other indicates the restoration of chiral symmetry. 

We should remind the reader that, in this paper, the probe branes only serve as the place 
where the string hangs. Our calculations are independent of the detailed brane configura- 
tions. 



3 Calculation of the momentum diffusion coefficient 
3.1 Preliminaries 

In [91 [10] the authors obtain the transverse momentum diffusion coefficient in the following 
way. Firstly, by analogy to classical theory for Brownian motion, they propose that 

«r = I [ dt (F(t)F(0) + F(t)F(0)), (5) 
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Figure 1: The dominant configurations of the D8 and anti-D8 probe branes in Sakai-Sugimoto 
model in deconfinement phase, (la) indicates the chiral symmetry is broken, and (lb) 
signifies the chiral symmetry is restored. If D8 and D8 are separated by a distance L at 
infinity, for T > 0.154/L, (lb) is the dominant phase, for < T < 0.154/ L, (la) is the 
dominant phase 

where F(t) is the transverse stochastic force acting on the probe quark. In M = 4 SYM, it 
takes the form 

F (t) = j d?xQ ] (t,x)T'Q(t J x)E'', (6) 

where E a is the field strength supplied by vector gauge fields and six scalar fields. We can 
define the Wightman correlation function and the Feynman correlation function about F(t), 
they are 

G(t) = ±(F(t)F(0)+F(t)F(0)), (7) 
G F (t) = -i(TF(t)F(0)). (8) 
And there is a relation between them in frequency space, 

G(u) = -lmG F (u). (9) 

So, Eq. flS} can be changed into 

kt = lim G(u) = — lim IhiGf(cg>). (10) 

uj^O ui— >0 

The ( ) denotes an average in the states which are composed of SYM and a moving quark. 
Using Wigner distribution function in QCD kinetic theory, the generating functional of the 
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above correlation functions can be written as the VEV of a Wilson loop. This loop is a closed 
contour in the complex time plane, specially, its x component should satisfy x = x +v(t c —t CQ ) 
and y component is equal to Syi(t c ) when t c lies on the real time axis, Sy^it^) when t c lies 
below the real time axis. This choice is determined by that the probe quark is traveling in 
x direction with velocity v, and 5yi(t), 5y2(t) are the fluctuations of quark's displacement in 
transverse direction acting as external sources coupling to the transverse stochastic force. 

The authors of paper J8[ [H [10] evaluate the VEV of the Wilson loop via AdS/CFT 
correspondence in ArfSs-Schwarzschild background, by finding out a string's classical action, 
requiring that the boundary of string's world-sheet is the Wilson loop. This is to say 

1 _ e iS NG [SyM = (W[S yi ,8y 2 }). (11) 



e iS NG [0,0] (W[0,l 

The closed time contour corresponds to the two boundaries of global AdSs-Schwarzschild 
space and the string stretches between the two boundaries. Finally, they obtain 

K T = v^AT 3 tt , (12) 

when v — > 0, 7 — > 1, Kt and the drag coefficient [2Tj [22] satisfy the Einstein relation. 



3.2 Calculation of the k t in S-S model 

In S-S model, E a appearing in the stochastic force term ([6]) should also include the con- 
tribution from K-K modes with mass of the order of QGP temperature. We will use the 
deconfinement phase background ([3]). For convenience, we use ut to scale dimensional coor- 
dinates and other parameter. 

, i¥ , , , n . fR\ 3/2 



/7/\3/2 / /? Y 

ds 2 =u\ (-J [-f{u)dt 2 + 5 ij dx l dx J +dxl] + ( 



u 2 dVt\+ 0,1 



/(«)J 



(13) 



where R = — — , f(u) = 1 — K, and now all the coordinates are dimensionless. Since this 
background is spherically symmetric , to define its Kruskal coordinates is a routine. They 
are 

U = e- 2u+ , V = e 2u ~, (14) 

where 

t — z*, (15) 
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and 



du 



(16) 



/W(|) 3/2 

But we only need the near horizon behavior of the Kruskal coordinates. In the near 
horizon limit, the metric becomes 

R 3 du 2 ' 



ds 2 ~ 3u 2 T R~ 3/2 



If we define 



1 - -)dt 2 + 

- u 9(i -a j 



^3/2 

P = — ~ — it, 2M 



+ i4#~ 3/2 [5 ij dx i dx j + dx\ + i? 3 ^] . (17) 
R 3/2 , . 



then this metric looks like Schwarzschild metric 



ds 2 ~ 3 M |i?- 3 / 2 



2M X , 2 dp 2 
-(1 )dt 2 + 1 



+u 2 T R- 3/2 [5 ij dx i dx j + dx\ + R 3 dQl] . (19) 



p ■ 1 - 2M/p_ 

The near horizon Kruskal coordinates are the same as in the Schwarzschild metric, 

U = -4Me-(^)/ 4¥ , V = 4Me( ( - r *'/ 4M , 

where 



(20) 



r* = p + 2Mln|p/2M- 1|. (21) 
In background (fTBl the string configuration with a constant velocity v in x 1 direction is as 



1/ ^ ^ 3/2 , 

x [t, u) = vt H z — In ■ 



\u - II 



vR 3 / 2 2u+l 
arctan ■ 



(22) 



3 " \J u 2 + it + 1 v 7 ^ ^3 
with other coordinates kept constant, where we have chosen the t and u coordinates to 
parameterize the string world-sheet. So the string configuration with a perturbation 5y(t, u) 
in the x 2 direction should be 



x M = (t, x (t, u), 5y(t, u),u, const). 
If we insert x^ into metric (fTBl . we get the induced metric on the world-sheet, 

3/2 



(23) 



d* 2 

U,J induced 



U \ ' 



R 



D4, 



' U T\ 7 ,2 



7 Z IT M 



^d^ 1 ~ ^ du2 , 1v 2 R%lurpdtdu 2 
H -3 h ddy(t,u) 



:i-S) 2 



it (1 



(24) 
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where we have restored the dimension of the coordinates, and 7 is the Lorentz factor. This 
metric can be simplified by performing coordinate transformation, 



- t 
t = - + 



R 



3/2 

D4 



7 37^ /2 



In 



V 



\lv 2 +v + 1 

— \/3 arctan -^—= 7 2//3 In 



V3 



where p = — , o 



u — u, 
then the metric becomes 



+ 7 2 / 3 v / 3 arctan 



2g 



V3 



""^induced 



U 



R 



3/2 



2 3 (Es£)^du 2 



l£ 3 



(25) 



(26) 



(27) 



Now we define iy^^ur = t, u^^ut = u, 5y{i^u)^ 2 ^UT = Sy(t,u), R^ 2 ^ut = Rda, then 
the metric looks like the original one (fT3|) . 



ds 2 = ( 7 2/3 M T ) 2 ( J R/n)- 3 / 2 



-f(u)dt 2 + (?) 3 |S- + d6y(t,u) 2 



(28) 



u f(u) 

with ut — > 7 2//3 mt, R —> R. A compelling characteristic is that from the point of view of 
string world-sheet, the horizon shifts to u — 1, or u=7 2 / 3 . This new horizon is usually called 
the world-sheet horizon. The near horizon Kruskal coordinates can be defined by 

U = -AMe-^/^K V = AMe^/^ 



where 



A R 3/2 „ n „-v i? 3/2 

p = u, 2M = 

H 3 ' 3 

r* = p + 2Mln|p/2M- 1| 



(29) 

(30) 
(31) 



In Figure 2, in the R patch of world-sheet Kruskal plane, the action of the small fluctu- 
ations 8y(t,u) is derived from the Nambu-Goto action 

S = (l ll U J )2 J dtduJl- f^{u)5y 2 + f(u)t 3 6y' 2 , (32) 

where "dot" denotes df, "prime" denotes 8%. Because 5y is small, we can expand action 
around Sy = 0, and keep up to the second order of Sy. 

S = / dtdu{l - \r\m 2 + \f{u)j- 5y' 2 ). (33) 



V 




Figure 2: The left figure is the space-time Penrose diagram in Kruskal coordinates. The 
right figure is the string world-sheet Penrose diagram in world-sheet Kruskal coordinates. 
The regions with yellow color in two figures, represent the same zone in space-time, so are 
the blue regions. The hyperbolas in the space-time Penrose diagram correspond to the 
world-sheet horizon. 



Note that the infinite part of the action is subtracted since it appears in the numerator and 
denominator of Eq. ffTTT) . To solve fluctuation 5y , we define 



2tt 



(34) 



where we have chosen to normalize Y(u = oo,u>) = 1, because y(£o) is the Fourier trans- 
formation of the boundary value of Sy(t,u). The Euler-Lagrange equation of small string 
fluctuations can be written as 



3u 2 



u" 



-d*Y~ 



2 u*R 3 



1) 



:Y* = 0. 



(35) 



This equation is solved by 



Y{u,Cj) 



u- 



F(l>,u) 



(36) 



53/2 



corresponds to in-falling fluctuation 



where F(lj, u) is a regular function of u. (1 — 4%) 
in the world-sheet horizon u = 1. The complex conjugate of this expression is also a solution 
of the differential equation ( 1351) and corresponds to out-going fluctuation in the horizon. 

Now, we have obtained the solution in the R patch of world-sheet Kruskal plane repre- 
senting the right part of the R patch of the space-time Kruskal plane. The right and left 



parts in the R patch are separated by the curve u 
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2/3 



as in Figure 2. But in order 
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to use gauge/string duality to obtain the generating functional of the Feynman correlation 
function, we need to know the solution defined in the whole Kruskal plane of space-time. 
To this goal, we will extend the solution in the right half of the R patch into other parts of 
Kruskal plane one by one. Firstly, we should extend the solution into the whole R patch of 
space-time Kruskal plane. In other words, this amounts to extending the solution from R 
patch of world-sheet Kruskal plane to the R and F parts of world-sheet Kruskal plane. In 
terms of near horizon Kruskal coordinates, the in-falling and out-going solutions in R patch 
of world-sheet Kruskal plane behave as 

in-falling: e~ iGlt Y{u, u) ~ e" 4 ^ ln ^, (37) 

out-going: e-^Y*(u,u) ~ e 4 ^ 1 *-^. (38) 

Because the Kruskal coordinates are global, actually, in the F patch of world-sheet Kruskal 
plane, near world-sheet horizon, we can also write down these two solutions as 

in-falling: e'^Y {u , u) ~ e^*™ lnV , (39) 

out-going: e~ iCli Y*(u,u) ~ e 4 ^ lnf> . (40) 

From above expressions, in the F patch of world-sheet Kruskal plane, we see that the in- 
falling solution is still effective, because V > 0, but U changes sign. So we will do an 
analytic extension for the out-going solution to make it an solution in the R and F patches 
of the world sheet Kruskal plane. Following Herzog and Son's prescription [23], the out- 
going solution should cross the horizon from the upper half of complex U, this results in 
that the out-going wave should picks up a factor e . Physically, this indicates that the 
out-going wave should be purely negative-frequency. We can repeat this process in the P 
and L patches. Having done this, the out-going solution picks up a factor e~ 47r£jM . Now, we 
have solutions defined in V > 0, V < parts of the world-sheet Kruskal plane. As we have 
demenstrated previously, the R and F patches of the world sheet Kruskal plane represent the 
R patch of the space-time Kruskal plane. The L and P patches of the world sheet Kruskal 
plane represent the L patch of the space-time Kruskal plane. So we know how these solutions 
behave in the R and L patches of the space-time Kruskal plane respectively. Near space-time 



10 



horizon, in terms of the space-time near horizon Kruskal coordinates, in the R patch of the 
space-time 

in-falling: e~^Y(u, u) ~ e - iitjMlnV , (41) 
out-going: e 4 ™™ e^Y* (u , u) ~ e ^ e - totflnV . (42) 
in the L patch of the space-time 

in-falling: e - iU Y(cu, u) ~ e - 4i " Mln( - y) , (43) 
out-going: e -^* e -^ty*(^^) ^ e -4^* e -4fo,Min(-v3_ ( 44 ) 

These expressions tell us that in the point of view of space-time, the in-falling and out-going 
solutions in the world-sheet are both in-falling waves. Moreover, the R patch solutions can 
be interpreted as the solution in both the R and F patches. The L patch solutions can be 
interpreted as the solution in both the L and P patches, because V > in both the R and F 
patches; V < in both the L and P patches. So far, we have obtained the solution defined 
in the whole V > part of Kruskal plane, and the solution defined in the whole V < part 
of Kruskal plane. They can be deduced from the following four different solutions defined in 
R and L patches. 

. e- M Y(u,u) in R j in R 

vr, i = \ yL,i=< .... / , (45) 

in L I e- luJt Y(oo,u) in L 

e- iajt Y*(u,u) in R ( in R 

VR, o = \ Vl,o=< .... , s • (46) 

inL e~ tuJt Y* uj,u) in L 



Following the Herzog and Son prescription [23], we look for linear combinations of these 
expressions that, close to the horizon, are analytic in the lower half of the complex V plane. 
Physically, this means that the in-falling wave should be purely positive-frequency. With 
this requiriement, the two linearly independent combination are: 

Vo = Vr, o + a VL, o, Vi = Vr, i + cxiVh, i ■ (47) 

where the a Q and can be determined from the near horizon behaviors of these solutions 
Eqs. flUD, (S2D, (|33D and (g3D 

a D = e ^M e -4^M ) (48) 
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These two solutions are used as basis for the linearized string fluctuations defined over the 
full (AdS) Kruskal plane 

duo 



y{t,u) 



2tt 



(a(u)y a (uo) + b(u)yi(u)) 



(50) 



The coefficients a(cu), b{Co) can be determined by the boundary values of the solutions. 
Because we have 

duo 



y(t,u = oo) 
y(i,u = oo) 



R 



duo 



2tt 



(51) 
(52) 



We obtain 



a(uo)=n(- yi (uo) + e 4 ™ M y 2 (uo)), (53) 

b(uo) = n{e^Uuj) - e^ M y 2 m- (54) 

where n = l/(e SlTwM — 1). Now we compute the boundary action in terms of the string 
solution: In (t, u) coordinates 



S 



B 



2ira' 



duo „ , . , v? 



R 



duo ii 3 



-f(u)—y(-uJ,u)duy(ud,u)- / —f(u) 1 -y(-uo,u)duy(uo,u) 



• (55) 



2tt R3' 

Notice that Co = ju, yi(uo) = 7~ 1 y 1 (<^), y 2 {Co) = 'y~ 1 y 2 (uo) } and using Eqs. (USD, (@SD, d3HD, 
( 149]) . (1531 . fl54|) and (J50~]l . this action can be expressed as 

_ ( 7 2/ V) 2 /"dw 
O 



27TQ; / 7 ,/ 27T 

{{n + l)f{u) T -Y*(-uo,u)duY(uo,u) - nf(u) 1 ~Y(-uo,u)duY*(ud,u) 
R 6 R 

~ 3 "3 

+ yi {-uo)y 2 {uo)e™l 2 h (-/(<*) "y*(-£, u)d & Y{u, u) + f(u)^-Y(-uo, u)d^Y*(uo, u) 

R 6 R s 

"3 "3 

+ y 2 (-a;)y2(w) (nf(u)-r-Y*(-uo, u)duY(cd, u)-(n+ l)f(u)—Y(-uo, u)duY*(uo, u) 

R A R A 



I U— >+OQ 



(56) 

From this expression we can read off the Feynman correlation by taking derivatives with 
respect to yi(-uo), yi{u). 

1 



= 9 ,'2/3 [(n + l)/(«)^W^«)-w/(u)-gry9fty*(^6) 
zna'^ z / 6 UT R 3 R 3 



(57) 
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where we have restored physical dimension of Gp- When w -> we can expand Y(uj,u) in 
a power series in Co and solve order by order 



Y(u,u) = (1 - -^y 2l£l " \l - 2iMcu(3ln(u/Vu 2 + u+l) 



U 3 ' 



-V3 arctan(2M/ v / 3 + l/VS) + ^— ) + 0(lu 2 ) 



(5f 



So 



3 7 l/3^2 

k t = - lim ImGV(w) = - 3 , (59) 



where we have used R 3 = 7 2 i? 3 . Using the parameter relation between string theory and 
gauge theory, we obtain 

16y^7T 7 1 / 3 AT 4 
«T = (60) 

where A is the the 't Hooft coupling of YM, gy M N c . When v — > if k,j< and the drag coefficient 
which is often denoted as t)d satisfy the Einstein relation, then the drag coefficient in S-S 
model should be 

Vn = ^(l^) 372 - (61) 
This is the result appearing in [19]. 



4 An estimation of q with respect to gluon energy loss 
in S-S model 

4.1 The initial energy- momentum of gluon in S-S model 

In |llj . S. Gubser and his collaborators proposed a simple but interesting idea to estimate 
the jet quenching parameter q, by computing how far an off-shell gluon propagates in the 
finite temperature Af = 4 SYM before it loses all energy and resolves into the medium. Their 
estimation is mainly based on an extension of BDMPS formalism [2U [251 126] from light-like 
parton to time-like parton, which is to replace the light cone distance L~, by \^2Ax, where 
Ax is the parton's in- medium space distance, often called penetration length. In short, this 
is 

1 L~ 2 1 

AE = -a s C R q— — > AE = -a s C R qAx 2 , (62) 
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and 



AAE 



In above expressions, a s is the strong coupling constant, and Cr is the color group SU(N) 
Casimir C^-R) evaluated in the parton's representation, for gluon Cr=N. Although this is 
not an exact calculation about q, it has a striking virtue appealing to us, that is their result 
q ~ 21GeV 2 /fm lies within the 3a range of averaged q [27], while other people's result [6] 
fails. The 3a range of averaged q is 



7 <^1 < «> < ^ 



(64) 



r r^j \j./ r^j r ' 

tm fm 

with lowest \ 2 & t (q) ~ 13GeV 2 /fm. S-S model is also a holographic dual to QCD, so we 
can ask whether their method still has this advantage in S-S model. 

To answer this question, we do the following estimate in S-S model. Firstly, we will 
introduce the following background for the simplicity of calculation, 

R 3 dy 2 . 



ds 2 = u 2 T {RyY^ [-(1 - y s )dt 2 + dx 2 + 



(65) 



y(i -y 3 ) 

This is obtained from the metric ( |T3i) . by replacing u with ^ and letting £4, ^4 be constant. 
Following S. Gubser's approach, a gluon is represented by a doubled string which rises from 
the horizon up to a minimum y^y, and then falls back down to the horizon as in Figure 3. 



Y=yuv 




Figure 3: The blue line shows the shape of a doubled string which represents a gluon. 
The pink line and the orange line represent the world-sheet light signal's trajectory and the 
massless particle's trajectory we will introduce respectively. We will use these two trajectories 
to give an estimation of the range of gluon penetration length. 
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Required to be stable, the gluon's initial state is constructed from the trailing string 
solution, the trailing string is a string moving with constant velocity in x 1 direction in the 
background fl65|) . and its shape is 

ln(l-y) Hl-(3y) ln(l - ay) 



x\t) 



vt + vRz 



+ 



+ 



3 

vR 2 



(l-a)(l-/3) (a-l)(a-(3) (J3 - a)(J3 - 1) 
ln(l - (3y) ln(l - ay) 



+ 



(66) 



a — (3 a — (3 

This solution is equal to (T22J), a = e 2m l 3 and /3 is a's complex conjugate, we prefer this 
expression because of its convenience for following calculation. So, if we insist to use world- 
sheet coordinates a a = (t,y), initially, the induced metric on the world-sheet is 

-h + v 2 -v 2 yR 3 l 2 /h 
-v 2 yR 3 ' 2 /h R 3 {h + v 2 - hv 2 )/h 2 y i 

where h = 1 — y 3 . The world-sheet current density of energy-momentum 
is 



g a/3 = u 2 T (Ry) 



(67) 



1 







(68) 



-h — v 2 + hv 2 v 
7ia'h{l-v2) y-hv 2 y 2 R- 3 ' 2 hvy 2 R~ 3 l 2 

where m=(0, 1, 2, 3) corresponds to (t,x). Usually, the related doubled string's energy- 



momentum in these four dimensions is 



Pn 



Ut 



— / dy^^P l n 



ut 



na' y/l - v 2 



dy 

hy 2 



7-^ I -h - v 2 + hv 2 v 



(69) 



Because dt, d% is the Killing vectors, so p m can be identified with the four- momentum of the 
gluon in the boundary gauge theory. But we are also confronted with the problem appearing 
in JV=4 SYM: the energy- momentum has a logarithmic divergence at y—1. Gubser gives an 
explanation for the appearance of this kind of divergence: this divergence is due to the fact 
that, to form the shape of a trailing string needs infinitely long time, and during this period, 
infinite energy-momentum has been transferred from the string to the medium, but it is still 
contained in the right hand side of Eq. f|69|) . So once this shape has been formed, the rest 
energy-momentum of the string should not include the energy-momentum transferred into 
the horizon. To make this subtraction, we compute 



fixedx 1 



- [ dy^-g[PL-(^- 

u t y WTV dy 



9[P m ( q l^i^m. 

_X 1 (-1,1,, 0,0 , 

na' y uv 



(70) 
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where the subscript x 1 indicates that the above integral is carried out along the x 1= conctant 
contour in the (t,y) plane. With Green theorem and P°. a =0, it is not hard to prove the 
difference between (1691) and flTOj) is the amount of energy-momentum we want to subtract 
from ( l69l . Now, we can show to the reader that this gluon is a time-like one. This is obvious, 
since 



4.2 Estimation of gluon penetration length 

Since this doubled string's tip does not attach to the boundary brane, this string will fall 
down toward horizon. At some moment, the tip will touch the horizon, from the beginning 
to this moment, the tip travels Ax in x 1 direction. Generally, Ax is function of y\jy and 
velocity v or 7. For a fixed energy gluon, Ax is function of yuv or v, for example, we choose 
Ax=Ax(v), then the maximum value of Ax(v) with respect to v can be interpreted as the 
penetration length of the gluon. Although, we know the initial shape of the doubled string, 
it is difficult to compute Ax in terms of the EOM of string, not to mention Ax max . But with 
the methods proposed by Gubser, we can find out a lower and an upper bound for Ax max . 

Firstly, we consider the lower bound. The initial shape of the doubled string can be 
interpreted as cutting the trailing string which attaches to the boundary brane, at y — ?/uv 
at some time. Meanwhile, a light signal is emitted from the cut. Since the disturbance 
arising from cutting the string cannot propagate more quickly than light, the string will 
keep its shape where the light has not even arrived, as if we did not make such a cut. 
The displacement of light in x 1 direction, denoted as Ax\ ow , should be the lower bound of 
Ax, then Axi ow , max serves as the lower bound of Ax max . The trajectory of the light signal 
can be determined from the light-like tangent vector of world-sheet metric, which satisfies 
as t increases, y increases. There is only one light-like tangent vector field meeting this 
requirement. It is 




l)] 2 > 0. 



(71) 




(72) 



So the light signal's 



trajectory is determined by 



dt 
dy 




v 2 y yl — f 2 



(73) 



h(h — v 2 ) y/y(h — v 2 ) 
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Solving this differential equation, we obtain 

2 



3 2 
t = i?2^3 



(1 -a) (1-/3) 



ln(l + 7 1/3 W) + 



( a _l)( a _/3) - 1-/3 7 V3^ 



In 



3 2 



In 



1 + a 7 1/3 v^ 



(/3-a)(/3-l) l-a 7 1/3 v^ 
ln(l-y) | ln(l - (3y) | ln(l -o/y) 



-J2* 



_(l-a)(l-/3) (a-l)(ar-/3) (/3 - a)(/3 - 1)_ 
ln(l-/3y) ln(l-ay) 



a — (3 



a — (3 



i?2 7 3 



3 2 

+-R273 



ln(l - 7 2/3 %) ln(l - 7 2/3 a?/) " 
(a - l)(a-/3) + (l3-a)(/3- 1). 
ln(l - 7 2 / 3 /3y) ln(l - 7 2 / 3 cn/) 



a — /3 a — (3 

where the a and (3 are the same as Eq. (I6"B"I) . Plugging Eq. ([71]) into Eq. 



(74) 

we obtain the 



orbit 



x\z) = u T ^R§ 7 § [ln(l + 7 1 / 3 V ^) + a ln(l + /3 7 1/3 + /3 ln(l + a^ 3 ^)]. 



So 



Aa;i ow =ae 1 (l) -x 1 (y uv ), 
where we restore the physical dimension of x l . For convenience, we define 



E 



utR 
na'E 

Ut 



3 ' 

2 



(75) 



(76) 



(77) 



(78) 



Now we should find the maximum value of Axi ow with a fixed E = 7 (— 1). If we define 



i = i 2/3 yuv 



(79) 



then Ax low is a function of £. Usually, to find the maximum value of Ax(£)i ow , we will first 
find a satisfying 

9 e Ax(0iow, e=C , = 0. (80) 

But this is only the point making Ax(C,)\ ow a local maximum or minimum and may not be 
the global maximum or minimum. In fact, there is only one satisfying Eq. flHUj) . and 
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A£(£*)iow is the global maximum, this is supported by numerical result. When E ^> 1, we 
find that can be expanded in terms of E^, it is 



= 0.38036 + 0. 207807£~5 + 0.07596-ETs - 0.781227E- 1 + 0(E 



-6/5> 



And 



Axi ow , max = 0.84423E5 - 0.80997 + 0.04178E" + 0.71431-ET* + Q{E~ 



(81) 



(82) 



We exhibit the comparison of analytic result and numeric result in Figure (4). We find that 
when E ^> 1, the analytic result indeed matches the numeric result well. 



X 



*■ n-upperbound 
a-upperbound 
*: n-lowerbound 
V* a-lowerbound 




Figure 4: The triangles represent the numeric result of upper bound of penetration length, 
the stars represent the numeric result of the lower bound of penetration length, The two lines 
are the analytic results of the upper and lower bound of the penetration length for E 3> 1 
and their extrapolations to E not satisfying E ^> 1. 

Having found the lower bound of penetration length of the gluon, next, we shall look for 
the upper bound. To find the upper bound, we use the following picture. When the tip of 
the doubled string begins to fall toward the horizon, it happens that a light-like particle is 
projected from the tip with 5-dimension velocity proportional to the 5-dimension velocity of 
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the tip. In other words, the massless particle's trajectory is tangent to the tip's trajectory 
at this point. The massless particle moves along the geodesic and will fall farther in x 1 
direction than the tip, because the tip is also pulled by the rest of the string besides the 
gravity. So we can perceive the massless particle's displacement in x l before it falls into the 
horizon, as the upper-bound of Ax which we define before, here we denote this upper-bound 
by Ax U ppcr. Then the maximum value of Ax upper should be the upper-bound of Ax max , for 
fixed energy. To obtain Ax upper , we solve the EOM of the massless particle in the black hole 
background (|65l) with the initial condition required before. As we know, a massless particle 
can be described by the following action: 

1 f 1 dX^dX v 
S = ~ 2 J H-G^^l (83) 

where e is a Lagrange multiplier, G^ u is metric (I65i) . In the following, let's work in a gauge 
where t] — y and consider the trajectory of the form 

X° = X°(y) X l = X\y) X 2 = X 3 = 0. (84) 

Then 

S = jdyC £= 4(%r 3/2 ( _ Mx o ? + {x r f + ^ (g5) 

where prime denotes d/dy. Since the Lagrangian does not contain X°, X 1 explicitly, we 
immediately form two conserved momenta 

UT h(Ry)-^x ' MRy)' 3 ' 2 * 1 ' faR , 

Po = Pi = • (86) 

e e 

The equation of motion of e is a constraint: 

Ut 



7) 



y 2 VPo~ h Pi 

Because of our metric signature, the -po is energy, and should be positive, so po is negative, 
for X ' is positive, corresponding that the massless particle falls toward black hole, finally, we 
choose plus sign in (1ST)) for the trajectory. Then the shape of the trajectories is determined 
by 

^ = X 1 ' = CPl ( Ry ^ /2 = ^3/2 Pl/PO ,ggx 

dy u T VVV 1 ~ h Pi/Po ' 

where pi/po is due to the initial condition, since p\ and po are conserved quantities. Because 
the tip of an open string or a doubled string must move at the speed of light, at the moment 
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when the tip is formed, its 5- dimensional velocity should be proportional to l fJ, \ y =y uv , which 
is l a dX^/da a \ y=yvv , satisfying G^W = 0. For l a, s definition, refer to (172]) . Then we can 
derive pi/po 

p 1 /p Q = v Vl - v 2 - y uv 

v 2 Vuv ~ VI -v 2 (l - yfj V ) 
Now we calculate how far the massless particle propagates in the X 1 direction before falling 
into the horizon: 

A* uppcr = -R^u T f dy ^ (90) 

The maximum of Ax upper = ^f^p for fixed E is depicted in the Figure 4. For E ^> 1, using 
the same method as before, we find that 

& = 0.32141 + 0.13949iH + 0.03882iH + 0(iH) (91) 

Ax upperjmax = 0.99033Et - 0.79951 + 0.49157iH + 0(E~^). (92) 

At this moment, we should remind the reader that the analytic results do not match the 
numeric results as well as in the lower-bound case, the deviation may come from the approx- 
imate method we adopt to find out £*. 

So far, we have got two bounds of the gluon penetration length. Using these two bounds 
of penetration length to make a rough estimation of jet quenching parameter will be carried 
out in the discussion section. 



5 Discussion 

Following Gubser's description, one is able to extract the jet quenching parameter from 
momentum diffusion constant kt, by 

4b = ^, (93) 

v 

where "B" indicates Brownian motion, because we prefer to interpret it as part of the jet 
quenching parameter, this part is due to the Brownian motion effect. Since strongly coupled 
S-S model gauge field theory is different from QCD, there exists considerable uncertainty in 
how to translate the results calculated in S-S model into quantitative predictions in QCD. 
To characterize this uncertainty, we recall that in J\f=4 SYM, the optimum scheme is 

Ta/-= 4 = T QCD /3 1/4 = 280/3 1/4 g 2 YM N c = 5.5 . (94) 
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The factor 3 1 / 4 comes from the requirement that M=A SYM and QCD are compared at 
the same energy density. Similarly, we also require that the S-S model gauge theory and 
QCD are compared at the same energy density, and then choose the following scheme for 
S-S model: 

Tss = Tqcd/C = 280MeV/C T dQCD = UOMeV g 2 YM N c = 5.5 . (95) 
In above expressions, 

C = o.gm^TW^) 173 , (96) 

and in scheme (195]) . £ ~ 1.43. The explicit calculation of the parameter ( will be given 
in appendix. We see that ( is an increasing function of Tqqd- This is in accordance with 
our previous argument that as temperature increase, more K-K modes will become active, 
for their masses are in tower of mTd, with m = 1, 2 • • • . Since more degrees of freedom 
contribute to the energy density, ( must increase accordingly to keep the energy density of 
S-S model plasma equal with that of QCD plasma. 
At this moment, we can calculate q B as following, 

q B ~ 5.4GeV 2 /fm^- = 1.85GeV 2 /fm- — . (97) 

Q 6 v v 

For charm quark, m c =1.4GeV, typical p c =10GeV/c 

q B ~ 3.5GeV 2 /fm. (98) 

This value of q B is larger than Gubser's result [10] q B = 1.4GeV 2 /fm and Liu's result [7] 
q = 0.86GeV 2 /fm. 

Now let us use the gluon penetration length to estimate the jet quenching parameter q. 
Following (195]) . using BDMPS formalism and setting a s = 1/2, we find 

E = 1.115.EGeV Ax max (E) = 0.236Ax( J E) max fm (99) 

q = GeVVfm. (100) 

Ai(£) max 



Because Ai(£) lowcr max < Ax(E) max < Ax(E) 



upper max) 



10.350AE „., 2/ „ A 10.350A.B „_ ra/ , , 1A1 , 
— — GeV 2 /fm < q < — — GeV 2 /fm. (101) 

Ax(-£/) U pp er max Ax(-£/)i ower max 
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A representative range of energies for hard gluons in the QGP produced at RHIC, is 5GeV < 
E < 25GeV. When we make quantitative estimate of q, we assume the energy of the gluon 
is 25GeV. The reason is that the higher the energy is, the longer the penetration length 
will be, and the ratio of radiative energy loss to collision energy loss should be larger, then 
AE appearing in BDMPS formalism may be roughly interpreted as the whole energy of the 
gluon, in short, in the following estimation, AE = 25GeV. Inserting this value of AE into 
Eq. ffTOTl) . we find 

89GeV 2 /fm < q < 106GeV 2 /fm. (102) 

This result is much far away from the experimental result 7GeV 2 /fm < q < 28GeV 2 /fm. 
The penetration length in S-S model seems too short, which may be a consequence of the 
extra hadronic degrees of freedom. 

If we insist that scheme (1951) should be the suitable one for comparing S-S model plasma 
with QCD plasma, the good performance of Gubser's method in J\f—4 SYM does not take 
place in S-S model. So there may be some unknown physical reasons that make Gubser's 
method break down when applied to S-S model. It is also possible that the numeric range 
of q is similar to that of the experiment may be a coincidence, so we could not require it to 
be a feature belonging to all holographic QCD models. But even if this method can not be 
always used to estimate q, the relation between gluon's energy and penetration length may 
be still meaningful. However, we do not know how to relate it to experimental observables. 

6 Appendix 

In this section, we will exhibit a detailed calculation of the parameter £ appearing in Eq. 

We recall that one of the first finite-temperature predictions of gauge/string duality is 
that of the thermodynamic potential of dual gauge theory in the strong coupling regime. 
The entropy is given by Bekenstein-Hawking formula S = A/ AG, where A is the area of the 
horizon, G is the ten-dimensional Newton constant. To evaluate A, we cannot use string 
frame metric (113j) . but should use the metric in Einstein frame [28]. The metric in Einstein 
frame is obtained from multiplying the string metric ( 1T31) by \f~g^e~^ ) is the dilaton. The 
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result is 

9/8 / p„ \ 15 / 8 

.(103) 



) \-f{u)dt 2 + 5 lJ dx i dx j + dxl] + f — j 

it_D4 / V U J 



du 2 1 

u 2 dttj + 



/(«) 

There are some parameter relations to be used, besides those given in (j3J) 

9ym = ^ 2 9JsT d (104) 
G = 8A S 2 Z S 8 , (105) 

where T5_^ is the Hawking temperature of metric (Tl3|) and (I103p . is the critical tempera- 
ture, A is the 't Hooft coupling, g a and l s are string coupling and string length respectively. 
Using these relations, we obtain the entropy density of dual gauge theory by, 

ss-s = I = (^)VAA^. (106) 

Applying the following thermodynamic relations between entropy density s, pressure P and 
free energy density F, 

dP = -dF = sdT, (107) 
we find the pressure and energy density are 

Pss = ~Ts s „s, (108) 
o 

e s -s = Itsss- (109) 
o 

It is straightforward to compute the sound speed of S-S plasma 

vl = \ (HO) 

This value of sound speed implies that S-S plasma is not a kind of conformal fluid, for 
conformal fluid, v 2 must be 1/3. We recall that the energy density of QCD is 

e Q cn = y N c T QCD, (HI) 
so we can deduce the relation between T$-s and Tqcd by demanding that es-s = € qcd, 

Ts-s = TgcDlOMAX^iTQCD/Td) 1 / 3 }- 1 . (112) 
Then we can extract ( from above expression 

C = O.dUX^iTQCD/To) 1 / 3 . (113) 

This is the Eq. flSSD- 
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